Natural Gauge Fields in Propagating Polariton Superfluids 
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We show that the interplay between the structure anisotropy and the energy splitting between 
the TE and TM modes of a microcavity leads to the appearance of a gauge field for a propagating 
polariton condensate. The field texture can be tuned by rotating the sample and ranges contin- 
uously from a Rashba to a monopolar field. In the linear regime, the monopolar field leads to a 
remarkable focusing effect. In the interacting regime, the spin-orbit coupling induces a breakdown 
of superfluidity. The spatially homogeneous flows become unstable and dynamically evolve into 
textured ground states such as stripes and domain walls. 



The creation of synthetic gauge fields has been seri- 
ously addressed in the last years [TH3] . In condensed mat- 
ter systems, spin-orbit coupling (SOC) plays a central 
role in the spin-Hall effect [5j [6] , topological insulators 
[7J |H] and semiconductor-based spintronics [9] . The cre- 
ation of synthetic magnetic fields in atomic Bose-Einstein 
condensates (BECs) has allowed the formation of vortices 
without stirring pQ, while the fabrication of a synthetic 
gauge field of the Rashba type leads to the formation of 
spin domains [2 j. In semiconductor microcavities, both 
photons and strongly coupled exciton-polaritons experi- 
ence a spin-orbit coupling. It originates from the splitting 
between TE and TM modes corresponding to an effective 
magnetic field [lOj [TT] with a quadratic dependence over 
k and a 20 azimuthal dependence. An additional split- 
ting between linear polarizations at k = corresponds 
to a constant magnetic field associated with the crys- 
tallography axes [12 . While the TE-TM field leads to 
several fascinating effects, such as the optical spin-Hall 
effect [T3HT5] or the acceleration of magnetic monopole 
analogs [16] [17], it cannot be integrated in the Hamil- 
tonian as a minimal coupling. In other words, alone, 
it does not exhibit the properties of a true gauge field 
satisfying gauge transformations. The dispersion of the 
elementary excitations of a polariton BEC in the pres- 
ence of the TE-TM splitting is anisotropic [18 , but the 
dynamics of the condensate is completely stable. The 
physical reason is that the TE-TM field does not break 
the time-reversal symmetry [19], contrary to what hap- 
pens in the Rashba/Dresselhaus field [20] [21]. The ef- 
fect of trapping has also been discussed in the context 
of SOC atomic BECs: As a result of the parametric in- 
stability occurring at k = 0, the ground state of the sys- 
tem is found to demonstrate either stripes or vortices, 
rather than the Thomas- Fermi profile [22 . In contrast, 
the TE-TM field in polariton BECs allow for a homoge- 
neous (phase mixed) ground state. 

In this Letter, we show that the situation can be dra- 
matically changed when the polariton condensate is put 
in motion. Under typical experimental situations, the 
condensate can be loaded at a controlled wavevector ko 
by simply adjusting the angles between the laser beam 




FIG. 1: (color online). Schematics of the experimental config- 
uration necessary to produce a controllable polaritonic gauge 
field in a microcavity. The angle between the crystallographic 
axis (the effective magnetic field ft) and the horizontal is ip/2 
(cp). A laser with a vertical inclination \ an d making an an- 
gle of with the horizontal produces a propagating polariton 
beam with wave- vector ko = (ljo/c) sinx(cos#, sin0), where 
cjo is the laser frequency and c is the light speed. 

and the normal and horizontal axes of the cavity (see Fig. 
[TJ. At a critical velocity v* = ftk*/ra, for which the static 
ft and TE-TM fields exactly compensate, the dynamics 
of the low- wavevector excitations is governed by a natural 
gauge field of the Rashba kind. In the present discussion, 
ft = fl(cos (/?, sin ip) describes an effective in-plane mag- 
netic field associated with the crystallographic axis [TP] , 
as illustrated in Fig. [I] Its orientation in the (k x ,k y ) 
plane can be tuned by simply rotating the sample. By 
investigating the properties of the collective excitations 
of a condensate polarized along the static field, we show 
that the gauge field makes the homogeneous superfluid 
flow unstable i) against transverse excitations in a limited 
region below k* and ii) against both transverse and lon- 
gitudinal excitations for wave vectors above fc*, depend- 
ing on the condensate polarization. These two instability 
mechanisms trigger the formation of spin stripes and spin 
domain walls. A remarkable regime occurs when the po- 
larization of the condensate is perpendicular to the static 
field. In such a case, a slowly propagating flow is unsta- 
ble, whereas the superfluid motion is recovered above a 
critical velocity. This leads to an unexpected breakdown 
of superfluidity in the subsonic regime. 
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We consider the motion of a bosonic fluid governed by 
the spinor Gross-Pitaevskii equation 



ih 



dt 



H + ai^ 1 "^^ + a 2 *^ffx^, 



(1) 



where \£ = (^ + ,^_) T , ol\ and 0L2 are the intra- and 
inter- spin interactions, and o x is the anti-diagonal Pauli 
matrix. The single particle Hamiltonian, valid for both 
weakly coupled photons and polaritons in the parabolic 
approximation, is written as 

h 2 v 2 



2m 



f3(d y + id x 



\-p{d y -id x y 
h 2 v 2 

2m 



(2) 

where f3 = h 2 /(4m r ) is the strength of the TE-TM field, 
mg and m t represent the longitudinal and transverse po- 
lariton masses, m r = m£m t /(m t — me) is the relative 
mass and m = (me + m t )/2 is the averaged polariton 
mass. The Hamiltonian ([2| describes an ideal polariton 
BEC, in the limit of the very long life-time, decoupled 
from the thermal bath. Diagonalization leads to the fol- 
lowing single-particle spectrum: 



h 2 h 2 / O 2 
— ± \ f3 2 k* - fink 2 cos(2# -tp) + ^-, 
2m V 4 



(3) 



where k = yk 2 +k 2 and k = /c(cos#, sin#). Due 

to the anisotropy of the TE-TM field, the spectrum 
([3| encodes very interesting features, as summarized 
in Fig. |2j Let us take cp — 0, for a moment. 
For k x < fc c i, with k c \ = y^ra^/ft 2 , the energy 
is a monotonic function of k y (Fig. |2}d)); on the 
other hand, for k x > k c \ the lower branch e-(k y ) 
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bends and the parametric process k y 

[2mk%/ (2m r - m) - 4Qmm r / (h 2 (2m r + m))] J 
serving both energy and momentum is possible (see 
Fig. (2^)). As a consequence, a transverse instability 
of the condensate becomes possible. More interestingly, 
the lower and upper branches cross each other at the 
wavevector k* = ±^/Q/2f3 (see Fig. [2^1)). At this specific 
point, the TE-TM and the applied fields compensate and 
the dispersion relation is locally isotropic and linear near 
q = 0, e± — e(fc*) ~ h 2 q 2 /(2m) ± with q = k — k*. 

This corresponds to dispersion of a system in the presence 
of spin-orbit coupling of the Rashba type. For arbitrary 
values of (/?, such compensation takes place at the prop- 
agation wavevector k* = ±^Q/2f3 (cos(/?/2, sin^/2). 
Around this point (\q\ < fc*), the polariton Hamiltonian 
reads 



2m 



(q 2 



c cr ■ q + K y (er x q) 



(4) 



and 



where k x = 2m\J JJU / h 2 cos cp/2 
2m^/pTt/h 2 sinip/2. In opposition to TE-TM spin- 
orbit coupling in Eq. (|2|, the latter describes a true 
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FIG. 2: (color online) Polariton dispersion relation for <p = 
(relative to Q/2). a) Dispersion along the x-direction showing 
the two polariton branches e± crossing at the magic point 
k* = y/Q/2/3. The propagation wavenumber ko defines a 
circle limiting the validity of the Hamiltonian Q . The critical 
point k c i satisfies e-(k c i) = e+(0) = f2. Panels b), c) and d) 
depict the dispersion along the ^/-direction, respectively for 
k x = {0.4,0.97,1.1} /im. In c), a parametric process k y : 
—> ±k p occurs between a X-polarized mode and a diagonal 
(D) and an anti-diagonal (aD) modes. In d), we observe the 
signature of the gauge field, with the two branches touching 
each other at k y = 0. Panel e) depicts the field texture of 
the off-diagonal terms in |2| in the reciprocal space and panel 
f) zooms around the compensation wavevector k*. We have 
considered the following parameters: m — 4.5 x 10 _5 m e , Q = 
0.2 meV and (3 0.1 meV/xm 2 . 



gauge-field, minimally coupled Hamiltonian, where k x 
and K y represent the magnitude of the monopolar- 
and Rashba- field coupling constants. The gauge-field 
texture S(k) for the particular case of propagation along 
a horizontal magnetic field (9 = (p = 0) is illustrated 
in panel e) and zoomed around k* in panel f) of 
Fig. [2] Because of the specific angular dependence 
of the TE-TM field, very different textures of the 
gauge field can be experimentally realized by simply 
rotating the sample. This results in more general 
field (or spin) textures, related by the transformation 
S'(k) = K((p)S[K(-(p/2)k\, where K(<p) is the 2 x 2 
rotation matrix (see the supplemental information in [23] 
for illustrative cases). Hereinafter, we consider the case 
cp = for definiteness. To illustrate the effect of gauge 
field Q in real-life experiments, we consider photons 
(polaritons) propagating against a defect at the magic 
wave- number kg = /c*e x . A cylindrical defect of size 
a, described by the potential U(r) = Uo exp(— r 2 /a 2 ), 
populates the states q± = ±27r/a. The particles of the 
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FIG. 3: (color online) Effects of gauge field for propagating 
polaritons. a) A flow passing an obstacle without fields (f3 = 
Q = 0). b) Focusing effect in the presence of the monopolar 
gauge field of Fig. [5]f) at the magic point /c*. We have used 
the same field parameters of the Fig. [2] in both situations. 
Panels c) and d) depict the circular polarization degree after 
a laser pulse of duration 10 ps and spot 5 /xm, with a polariton 
life-time of 100 ps. c) Vertical incidence (x = 0), showing the 
spin pattern related to the spin eigen-states of the TE-TM 
field [13] . Panel d), polaritons propagating at ko = fc*e x 
(x — 10° for GaAs) exhibiting the two-fold spin pattern of in 
the presence of a gauge field. 



mode q y = q+ (q y = g_) are pulled down (pushed up) 
(see [23 for details). In fact, the defect focuses the 
beam, as if light was "bent" by it (see Figs. ^) and 
[SJd)). This is quite remarkable, since focusing is obtained 
with an impenetrable defect. Another signature of a 
gauge-field, which can also be very easily observed in 
basic experimental setups, is the conversion of a four-fold 
spin pattern, associated with spin Hall effect due to the 
TE-TM field [13 into a two-fold spin pattern. In Figs. 
[3]c) and[3]d), we depict the circular polarization degree 
p c = (n + — n_)/(n + — n_), with n± = \ip±\ 2 , f° r the 
case of a cw pulse of 10 ps pumping polaritons at k = 
and at the magic wave- vector ko = fc*e x . 

To include the interactions, we consider the dy- 
namics of a polariton condensate propagating in the 
vicinity of the "magic" point fc*. The collective 
excitations of a polariton BEC propagating with 
velocity vq = foko/m can be investigated by lin- 
earizing with the Bogoliubov prescription, \l/ = 



3 i(ko-r— fit/K) 



i(q-r— u;£) 



-i(q-r— out) 



*o + E q (u q e 

[23] , for a linearly polarized condensate, \l>o = 
(n + ,n_) T = no(e^,e _zr? ) T , with r] being the po- 
larization angle, that corresponds to a state with 
magnetic energy E = -(1/2) / dr^Qa x ^. Therefore, 
a condensate is said to be X- (Y-) polarized for 77 = 
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FIG. 4: (color online) Transverse (T) and longitudinal (L) 
instabilities for a X- polarized condensate with a blue-shift 
of am = 1 meV and cxi — — 0.2oa. The solid and dashed 
lines respectively represent the real and imaginary parts of e± . 
Panel a) k c i < ko < k*; panel b) ko > k*. Stability diagram 
for X- (panel c)) and Y-polarized (panel d)) condensates. The 
black solid line represents the velocity v* = hk*/m, and the 
red/light grey (blue/light grey) dashed lines correspond to 
Vcl = hk c i/m (v C 2 — hk C 2/m). 



(r] = 7r/2), and the associated pseudo-spin is So = noe x 
(So = — noe x ) [28]. The chemical potential contains the 
interaction and "magnetic" energy of the condensate 



/i = (a 1 + 0^2)^0 — ^ (p 1 ~ cos 2( 



I cos 277. (5) 



We consider a BEC propagating with wavevector kg = 
koe x and for which the non-linearity is large enough so 
the critical Landau velocity for the breakdown of super- 
fluidity is larger than = hko/m. Let us discuss the 
case of a X-polarized BEC first (X is the ground state at 
q = 0). For fco < fc c i, the dynamics of the condensate 
is stable, characterized by two modes e±, which, in the 
long-wavelength limit q^ < 1 read 



where eo = 60 (fco, ^, 



c v -q y , e + {q y ) a 



(4- 



x 2^2\ X / 2 



Qx 
V 2 2 



1/2 



(6) 



is the gap energy at the origin and 



c± y = c± y (ko, O, /?) represent speeds. For fco > fc c i, the 
two modes in the longitudinal direction become acous- 
tic, i.e. e± ~ c±q x . Interestingly, along the transverse 
direction the lower mode becomes flat with a positive 
imaginary part (fig [4] a)), announcing the onset of an in- 
stability at q y =0, as a consequence of the parametric 
process of Fig. [l]c). 

This purely transverse (T-) instability takes places in 
the interval < < k v . The most unstable mode is 
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FIG. 5: (color online) Snapshots of the diagonal polariza- 
tion degree pr> — (tld — n a D) / {ud + n a D) after the onset 
of instabilities in a X-polarized BEC. v/co — 0.45 (a) and 
v/co = 0.6 (b) . We have considered the following parame- 
ters: m = 4.5 x 10~ 5 m e , ^ 0.2 meV and f3 = 0.1 meV/im 2 , 
ain — 1 meV and «2 — — 0.2«i. 

achieved for \q y \ = ft p /2, at the point where the imagi- 
nary part of dispersion reaches its maximum (see Fig. [I] 
c)). Finally, for ko > k*, the dynamics becomes unstable 
along the longitudinal direction as well (T+L-instability) 
as shown the[4]b). This is a consequence of the paramet- 
ric process shown in Fig. [l]d). The most unstable mode, 
however, is now observed for (q Xl q y ) — (0,ft p /2), similar 
to the zero wave-vector instability discussed for atomic 
BECs (see e.g. Ref. [26 ). The substantial difference 
is that here this instability is only achieved for a prop- 
agating (and not for a static) condensate. The phase 
diagram summarizing these results is depicted in Fig. [4] 
c). Increasing the velocity, the system evolves from su- 
perfluid to transversally, transversally and longitudinally 
unstable to finally reach the supersonic regime. 

The case of a Y-polarized BEC is summarized in Fig. 
[4] d) , and is even more intriguing. For very low wave- 
vectors, namely ko < fc c i, the spin is anti-aligned with 
the static field (S = — noe x ) and the system is T+L- un- 
stable. Above the magic value /c*, the imaginary part 
of the energy vanishes in the x direction, and the sys- 
tem exhibits a purely transverse instability. At wave- 
vectors higher than the second critical wave- vector k2 C = 
^/2m t Q/ti 2 , the system becomes stable. This counter- 
intuitive result simply states that superfluidity breaks 
down when the flow velocity is decreased or when particle 
density is increased. 

To confirm the previous results, we performed numeri- 
cal simulation using the spinor Gross-Pitaevskii Eq. 0. 
For a X-polarized condensate, we observe the formation 
of stripes along the transverse direction for ko > fc c i, il- 
lustrating the onset of a transverse instability (Fig [5] a)) 
and the formation of spin domain walls for ko > k* by 
the T+L- instability mechanism (see Fig. [5]b)). Repeat- 
ing the calculations for a Y-polarized BEC, the former is 
observed for k* > ko > k C 2 while the latter is obtained 
for ko < k* (see the supplemental information [23 ). At 
late stages, the instability saturates and the domain wall 
break into pair of half- vortices, the elementary topologi- 



cal defects in spin-anisotropic BECs [16] (a movie is pro- 
vided in [23 ). For both polarizations, we observe that 
the formation of spin domain walls is much faster than 
the formation of stripes, which can be understood from 
the polarization of the excitations. The onset of trans- 
verse (longitudinal) instability is the consequence of a 
parametric process which populates the modes q y = ±k p 
(q x = ±K X ) from the initial state at (q x ,q y ) = (0,0). 
The efficiency of the process is proportional to cos 2 (2 At?) , 
where A77 is the relative polarization angle between the 
initial and final states [24] [25] . From Fig. [I] we observe 
that a process of the type q y : —> ±k p occurs between an 
X- and two almost diagonal (D and aD) states, for which 
At? ~ 7r/4 and cos 2 (2 At?) <C 1. The longitudinal instabil- 
ity develops faster, since the process q x : — >• ±n x con- 
serves polarization and, therefore, cos 2 (2 At?) = 1. Our 
results reveal a curious difference in respect with atomic 
BECs, where exotic phases are formed in static systems 
[22j [26j [27] . Here, instead, the gauge field is formed due 
to the balance between the TE-TM and the static fields, 
for which propagation is absolutely required. The config- 
uration we propose is not only realistic, but is realized in 
any microcavity where both static and TE-TM induced 
fields are always present. 

To conclude, we would like to discuss the possibility 
of engineering a quadratic spin-orbit term, similar to the 
TE-TM field, for the case of atomic BECs. The general 
scheme to synthesize a gauge field of the Rashba type 
is via Raman coupling between of the Zeeman splited 
states rriF = —I and mp = 0, belonging to a certain tar- 
get state F = 1, in such a way that the state mp = +1 is 
far detuned. The effective gauge field, linear in fc, results 
from a simple rotation of the Pauli matrices around to 
their axis (see Ref. [2] for details). If the Raman cou- 
pling targets the spin states mp = ±1, with a far-detuned 
state ttif = 0, the rotating wave approximation and the 
reduction to a two-level problem would imply a double 
rotation of the Pauli matrices, and therefore a spin-orbit 
coupling quadratic in k and similar to our hamiltonian 
([2| would be possible for atoms. The achivement of a 
'TE-TM" field for atomic condensates would be very in- 
teresting and would motivate the discussions between the 
two communities, leading to the observation of interest- 
ing topological effects like the nonlinear spin Hall effect 
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